Tunnel Magnetoresistance of a Single-Molecule Junction 



Alireza Saffarzadeh 1,2 EI 

1 Department of Physics, Payame Noor University, Nejatollahi Street, 159995-7613 Tehran, Iran 
2 Computational Physical Sciences Laboratory, Department of Nano-Science, 
Institute for Research in Fundamental Sciences (IPM), P.O. Box 19395-5531, Tehran, Iran 

(Dated: December 26, 2008) 

Based on the non-equilibrium Green's function (NEGF) technique and the Landauer-Biittiker 
theory, the possibility of a molecular spin-electronic device, which consists of a single C6o molecule 
attached to two ferromagnetic electrodes with finite cross sections, is investigated. By studying the 
coherent spin-dependent transport through the energy levels of the molecule, it is shown that the 
tunnel magnetoresistance (TMR) of the molecular junction depends on the applied voltages and the 
number of contact points between the device electrodes and the molecule. The TMR values more 
than 60% are obtained by adjusting the related parameters. 



I. INTRODUCTION 

Traditional magnetic tunnel junctions use inorganic in- 
sulators as spacers 0, 0. However, recent experimen- 
tal data have clearly shown that organic molecules can 
serve the same purpose and a rather large TMR can be 
found [1 H H [|, S 0, II- In these experiments, spin 
polarized transport through molecular layers sandwiched 
between two magnetic layers has been demonstrated for 
systems involving carbon nanotubes [H, other organic 7r 
conjugate systems 0, [H, H, 0], molecular bridges [8[, and 
self-assembled organic monolayers @ . Organic materials 
have relatively weak spin-orbit interaction and weak hy- 
perfine interaction, so that spin memory can be as long 
as a few seconds [H| . Such features make them ideal for 
spin-polarized electron injection and transport applica- 
tions in molecular spintronics. 

Spin-polarized transport through organic molecules 
sandwiched between two magnetic contacts has also been 
recently investigated theoretically [ill, E2, EH EH EE 
El, EH El, El mi M, SH H- Many of these calcula- 
tions, based on density functional theory or tight-binding 
model, have shown that by changing the magnetic align- 
ment of the contacts one can substantially affect the elec- 
tronic current in the molecular devices. For instance, in 
a single benzene- 1-4-dithiolate molecule sandwiched be- 
tween two Ni clusters, parallel magnetic alignment led to 
significantly higher current (about one order of magni- 
tude) than antiparallel alignment, suggesting the possi- 
bility of a molecular spin valve [111 ]. 

Among many types of molecules, the fullercnc Ceo 
which is one of the most well-investigated organic semi- 
conductors, is suitable as a molecular bridge in magnetic 
tunnel junctions, because its lowest unoccupied molec- 
ular orbital (LUMO) is situated at relatively lower en- 
ergies in comparison with the other organic molecules. 
In recent years, the electrical and magnetic properties of 
C60-C0 nanocomposites, where Co nanoparticles are dis- 
persed in Ceo molecules have been studied 0, |2f| [2(| , 



and the maximum TMR ratio of about 30% at low bias 
voltages was reported. Recently, He et al. [13], studied 
the spin-polarized transport in Ni/Cgo/Ni junction using 
density functional theory and the Landauer-Biittiker for- 
malism. They showed that the binding sites of Ni on Ceo 
molecule play a crucial role in the transport properties of 
the system and a large value for junction magnetoresis- 
tance was predicted. We believe, however, that no theo- 
retical study on spin-polarized transport through a single 
Ceo molecule and its TMR effect, based on tight-binding 
method, has so far been reported. 

In this paper, we use a single Cgo molecule in between 
two ferromagnetic (FM) electrodes, instead of the usual 
insulator layer as in the standard TMR setup, to inves- 
tigate the possibility of a single-molecule spintronics de- 
vice. An interesting feature of the Ceo molecule is the 
emergence of a quantum loop current which is related 
to the degeneracy of the energy levels of the molecule 
and its magnitude can be much larger than that of the 
source-drain current [28| . This feature is useful for spin- 
polarized transport through the molecule. Our approach 
indicates that the spin currents in such a molecular junc- 
tion are mainly controlled by the molecular field of the 
FM electrodes, the electronic structure of the molecule 
and its coupling to the electrodes. 



II. MODEL AND FORMALISM 

We calculate the spin currents through a single Ceo 
molecule, sandwiched between two semi-infinite FM elec- 
trodes with simple cubic structure and square cross sec- 
tion (x-y plane). The model of such a structure is shown 
schematically in Fig. 1(a) and 1(b). Since the electron 
conduction is mainly determined by the central part of 
the junction, the electronic structure of this part should 
be resolved in detail. It is therefore reasonable to decom- 
pose the total Hamiltonian of the system as 

H = H l + H c + H r + Ht ■ (1) 
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The Hamiltonian of the left (L) and right (R) elec- 
trodes is described within the single-band tight-binding 
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the absence of FM electrodes is expressed as 
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FIG. 1: Schematic view of the FM/Ceo/FM molecular junc- 
tion at zero voltages for (a) parallel, and (b) antiparallel align- 
ments of the magnetizations in the FM electrodes, (c)-(e) The 
three different ways of coupling between the Ceo molecule and 
the magnetic electrodes used in this work. Black circles show 
the position and the number of couplings. 



approximation and is written as 
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, T (cf„,o-) creates (destroys) an electron with 
spin a at site i in electrode a (=L, R), the hoping pa- 
rameter ti a t j a is equal to t for the nearest neighbors and 
zero otherwise. Here, eo is the spin independent on-site 
energy and will be set to 3t as a shift in energy, — a ■ h a is 
the internal exchange energy with h a denoting the molec- 
ular field at site i a , and a being the conventional Pauli 
spin operator. The Hamiltonian of the Cgo molecule in 
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where d\ c a (di c ,a) creates (destroys) an electron with 
spin a at site i of Ceo and Ci c is the on-site energy and 
will be set to zero except in the presence of gate voltage 
Vq that shifts the energy levels of the molecule and hence 
€i c = Vg ■ The hoping strength ti Jc in Ceo molecule de- 
pends on the C-C bond length; thus, we assume different 
hoping matrix elements: t\ for the single bonds and ti 
for the double bonds. Most recently, we have shown that 
the effect of bond dimerization may considerably affect 
the electron conduction through the molecule under suit- 
able conditions [29]. Finally, Ht describes the coupling 
between the FM leads and the molecule and takes the 
form 
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The hopping elements U a j c between the lead orbitals 
and the ir orbitals of the molecule are taken to be t'. In 
this study we assume that the electrons freely propagate 
and the only resistance arising from the contacts. This 
means that the transport is ballistic 30]; therefore, we 
set t' — Q.5t, because the value of t' should not be smaller 
than the order of t. On the other hand, we assume that 
the spin direction of the electron is conserved in the tun- 
neling process through the molecule. Therefore, there 
is no spin-flip scattering and the spin-dependent trans- 
port can be decoupled into two spin currents: one for 
spin-up and the other for spin-down. This assumption is 
well-justified since the spin diffusion length in organics is 
about 4 nm [10] and especially in carbon nanotubes is at 
least 130 nm Q], which are greater than the diameter of 
Ceo molecule (~ 0.7 nm). 

Since the total Hamiltonian does not contain inelas- 
tic scatterings, the spin currents for a constant bias volt- 
age, V a , are calculated by the Landauer-Biittiker formula 
based on the NEGF method [H: 



/OO 
T CT (e ! y a )[/(e- Mi )-/(e- Mfl )]de 



(5) 



where / is the Fermi distribution function, [Il.r = 
Ep ± \eV a are the chemical potentials of the electrodes, 

and T a (e,V a ) = Tr[r L ^G a r R ^Gl.] is the spin-, energy- 
and voltage-dependent transmission function. The spin- 
dependent Green's function of the Ceo molecule coupled 
to the two FM electrodes (source and drain) in the pres- 
ence of the bias voltage is given as 

G a {e,V a ) = [el-H c -t Lta {€-eV a /2)-t R ^(e+eV a /2)]- 1 , 

(6) 

where £j, ](T and T, R (J describe the self-energy matrices 
which contain the information of the electronic structure 
of the FM electrodes and their coupling to the molecule. 
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These can be expressed as E a)0 -(e) = Tc, a ga,a(t)T a ,c 
where f is the hopping matrix that couples the molecule 
to the leads and is determined by the geometry of the 
molecule-lead bond. g a _ a are the surface Green's func- 
tions of the uncoupled leads i.e., the left and right semi- 
infinite magnetic electrodes, and their matrix elements 
are given by 
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a ■ h Q + e-(k) 
where r; = (ac», 3/i, k = (l x ,l y , k z ), z = e + i8 
2V2 
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Here, Z X)1 , (= l,...,N x , y ) are integers, k z e [— J, f ], 
and iVg with (3 = x, y, z is the number of lattice sites 
in the (3 direction. Note that N x and N y correspond 
to the number of atoms at the cross-section of the FM 
electrodes. Using E a)tr , the coupling matrices r aj(T , also 
known as the broadening functions, can be expressed as 
T a . a = -2Im(E Qj(J ). 

In the semi-infinite FM electrodes described by the 
single-band tight-binding model, only the central site at 
the cross section is connected to the molecule. Further- 
more, when the molecule is brought close to an elec- 
trode, the bonding between them will depend on the 
molecule orientation. This orientation can be such that 
only one carbon atom, a pentagon or a hexagon of the Cgo 
molecule be in contact with the leads [see Figs. l(c)-l(e)]. 
Therefore, one can expect different conduction through 
the molecule, which arises due to the resonant tunneling 
and the quantum interference effects. Our approach, as a 
real-space method, makes it possible to model arbitrarily 
the number of contacts. In this regard, the core of the 
problem lies in the calculation of the spin-dependent self- 
energies El j(T and E_R jCr . In the case of contact through 
a single carbon atom of the molecule, only one element 
of the self-energy matrices is non-zero. However, for the 
transport through opposite pentagons or hexagons, 25 
or 36 elements of the self-energy matrices are non-zero, 
respectively. 

The total charge current is given by I = Jf + Zj, . In 
this case, we calculate the TMR ratio from the usual 
definition: TMR = (I p — I a )/I p , where I p ^ a are the to- 
tal currents in the parallel and antiparallel alignments of 
magnetizations in the FM electrodes, respectively. 



III. RESULTS AND DISCUSSION 

We now use the method described above to study 
the coherent spin-dependent transport and magnetore- 
sistance effect of FM/Cgo/FM molecular junction. We 
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FIG. 2: The total currents at Va=0.0 V for the parallel (a) 
and antiparallel (b) alignments of magnetizations in the FM 
electrodes, and (c) the TMR as a function of applied voltage. 



have done the numerical calculations for the case that 
the direction of magnetization in the left FM electrode is 
fixed in the +y direction, while the magnetization in the 
right electrode is free to be flipped into either the +y or -y 
direction. We set |h Q |=1.5 eV, t=l eV, t\ — t, t 2 = 1.1 1, 
N x = N y = 5, and T=300 K in the calculations. 

In Fig. 2 we show the current-voltage (I-V) charac- 
teristics (in the parallel and antiparallel alignments) and 
also the TMR as a function of applied voltage for three 
different ways of coupling between the Ceo molecule and 
the FM electrodes. These three cases were chosen as the 
most probable experimental orientations s 3lj . As can be 
seen, the I-V curves show a steplike behavior which in- 
dicates that a new channel is opened. The values of cur- 
rents in the parallel configuration are nearly two times 
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larger than that in the case of antiparallel one. The 
difference between I p and I a is related to the asymme- 
try of surface density of states (SDOS) of the FM elec- 
trodes for spin- up and spin-down electrons [see Figs. 1(a) 
and 1(b)] [33], and the quantum tunneling phenomenon 
through the molecule. Here, we have assumed that only 
itinerant sp-like electrons contribute to the tunneling cur- 
rent, which is reasonable when the distance between the 
two FM electrodes is greater than 0.5 nm [33| . 

In the parallel alignment, minority electrons go into 
the minority states by tunneling through the molecule 
and there is no asymmetry in the SDOS. If, however, the 
two FM electrodes are magnetized in opposite directions, 
the minority (majority) electrons from the left electrode 
seek empty majority (minority) states in the right elec- 
trode, that is, asymmetry in the SDOS. In fact, in the 
antiparallel alignment, the tunneling currents for both 
spin channels are asymmetric with respect to voltage in- 
version (not shown). Consequently, the parallel arrange- 
ment gives much higher total current through the Ceo 
molecule than does the antiparallel arrangement. This 
difference in the total currents is the origin of TMR ef- 
fect which has been shown in Fig. 2(c). The TMR ra- 
tio has it maximum value (more than 60%) at low bias 
voltages. With increasing the applied voltage, we first 
observe that the TMR decreases. Such a behavior is sim- 
ilar to the conventional magnetic tunnel junctions [34| . 
With further increase in the bias voltage, the TMR ra- 
tio increases and then, the reduction and enhancement 
of TMR are repeated. Recently, using Cgo-Co nanocom- 
posites, Miwa et al. [26| measured a magnetoresistance 
value as large as 18% due to the spin-polarized tunneling 
of carriers between Co nanoparticles via Cgo molecules. 
Our bias-voltage dependence of TMR is qualitatively in 
agreement with their results. A similar bias-dependent 
behavior in which TMR ratio increases with applied volt- 
age, using a EuS spin-filter tunnel barrier in conventional 
junctions, has also been reported (35[. 

In order to investigate the other features of the junc- 
tion, we show in Fig. 3 the effects of gate voltage on 
the total currents and the TMR in the cases of single 
and multiple contacts. When the gate voltage is zero, 
we showed in Fig. 2(c) that the maximum values for the 
TMR are obtained at low applied voltages (V a < 0.1 
V). In such cases, the number of states contained in 
the energy window between /i£ and /ir is zero, because 
the energy window lies in the highest occupied molec- 
ular orbital (HOMO)-LUMO gap of the C 60 molecule 
[Figs. 1(a) and 1(b)], where there is no molecular level. 
Therefore, the current flow mechanism is tunneling. For 
this reason, the voltage dependence of the TMR effect at 
low voltages is similar to that of the conventional mag- 
netic tunnel junctions. It is worth mentioning that in the 
selected voltage interval, we did not observe significant 
changes in the transmission spectra with increasing the 
bias voltage. 

Applying a gate voltage shifts the molecular levels rela- 
tive to the Fermi level of electrodes, and hence the trans- 
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FIG. 3: The total currents at 14=0.1 V for the parallel (a) 
and antiparallel (b) alignments of magnetizations in the FM 
electrodes, and (c) the TMR as a function of gate voltage. 



mission coefficients may significantly vary [291 ] . Since the 
currents depend on the molecular density of states lying 
between and {j,r, then, by increasing the gate voltage, 
HOMO or LUMO peak moves inside the energy window 
and the total currents I p and I a increase. In this case, for 
either negative or positive gate voltages, the current flow 
mechanism is resonant tunneling, and a peak appears in 
the current curves. Figures 3(a) and 3(b) show that the 
currents vanish in a wide range of gate voltages, because 
in these voltages there is no resonant level inside the en- 
ergy window. In both figures, the peaks appear at the 
same voltages which confirm that the appearance of cur- 
rent peaks is due to the molecular levels. However, the 
height of peaks is different which is due to the difference 
in the spin-dependent SDOS of the FM electrodes. The 
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results suggest that the Ceo molecule is an interesting 
candidate for operation of devices as a nano-scale cur- 
rent switch. Also, with increasing the number of contact 
points between the device electrodes and the molecule, 
the interference effects around these points become im- 
portant, some resonances might completely disappear, 
and the spin current changes. This is the reason of dif- 
ference between the currents in the single and multiple 
contacts [see Figs. 2 and 3] [H, M H3 ■ 

IV. CONCLUSION 

Using the NEGF method and the Landauer-Buttikcr 
theory, we have investigated the possibility of making 
a Ceo-based magnetic tunnel junction. We have shown 
that coupling between the single Ceo molecule and mag- 
netic electrodes in FM/Cgo/FM structure, produces large 
magnetoresistance effects (greater than 60%) that can be 
modulated by using the molecule orientation, bias and 



gate voltages to control the spin-polarized transport. The 
present study advances the fundamental understanding 
of spin-dependent transport in molecular junctions and 
suggests that the Ceo molecule is an interesting candi- 
date for application in the magnetic memory cells and 
spintronic devices. 

Throughout this study, we have ignored the effects of 
inelastic scattering and the magnetic anisotropy of the 
FM electrodes due to the reduced dimensionality and the 
geometry of the system. These factors can affect the 
spin-dependent transport and hence, another improved 
approach is needed for more accurate results. 
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